NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


TECHNICAL NOTE 2249 

THE SPANWISE DISTRIBUTION OF LIFT FOR MINIMUM 

INDUCED DRAG OF WINGS HAVING A GIVEN LIFT 

AND A GIVEN BENDING MOMENT 

By Robert T. Jones 

Ames Aeronautical Laboratory 
Moffett Field, Calif. 







NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 22^9 

THE SPANWISE DISTRIBUTION OF LIFT FOR MINIMUM INDUCED 
DRAG OF WINGS HAVING A GIVEN LIFT AND A 
GIVEN BENDING MOMENT 
By Robert T. Jon.es 


SUMMARY 


The problem of the minimum induced drag of wings having a given lift 
and a given span is extended to include cases in which the bending moment 
to be supported by the wing is also given. As in the classical problem 
of induced drag, the theory is limited to lifting surfaces traveling at 
subsonic speeds. It Is found that the required shape of the downwash dis- 
tribution can be obtained In an elementary way which is applicable to a 
variety of such problems. Expressions for the minimum drag and the cor- 
responding spanwise load distributions are also given for the case in 
which the lift and the bending moment about the wing root are fixed while 
the s pan Is allowed to vary. The results show a 15— percent reduction of 
the induced drag with a 15 — percent increase in span as compared with 
results for an elliptically loaded wing having the same total lift and 
bending moment. 


INTRODUCTION 


In the problem of minimum induced drag as originally treated by Muhk. 
(references 1 and 2) the span of the wing and the total lift, were supposed 
to be given anfl the distribution of lift over the span resulting in a min- 
imum of drag was sought. The solution of this problem thus provided a 
convenient lower bound for the induced drag of a wing of given dimensions. 

In the practical design of wings the requirements for low induced 
drag and the requirements for structural strength are opposed. Here the 
bending moment developed by the lift becomes an important considera- 
tion — more important in many cases than the actual spanwise dimension 
of the wing. Such considerations lead to the problem of determining the 
minimum drag with limitations imposed on the bending moment as well as 
on the total lift. It is the purpose of the present paper to show how 
the methods of the earlier analysis can be extended in a very simple way 
to the solution of problems involving the bending moment of the load 
distribution. 
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A complete list of symbols employed in the analysis will be found in 
the appendix. 


GENERAL FORMULAS FOR LIFT, DRAG, AND BENDING MOMENT 


Reference may he made to the original papers of Prandtl and M uhk 
(references 1 and 2), or to any of the standard text books on aerodynamics, 
for the fundamental developments of wing theory which form the basis for 
the calculations of induced drag. In these developments the over-all lift 
is given by 

n +s 
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( 1 ) 


and the drag is given by 
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In these formulas the wing span is supposed to extend along the y axis 
between — s and +s, T is the local circulation or vortex strength, and 
Y is the constant velocity of flight. The induced downwash velocity w-^ 
is variable along the span and is connected with the vortex distribution • 

T (y) through the relation 



( 3 ) 


With this value for wq the expression for the drag may be converted to 
a double integral involving the spanwise distribution of lift as repre- 
sented by the circulation strength T 




r(y)r , (r]) 

y-*i 


dy drj 


( 4 ) 


This integral may be reduced to a more symmetric form if it is integrated 
by parts on the supposition that T falls to zero at the wing tips. Thus 1 



n + 3 r(y)r(q) 


dy dq 


( 5 ) 


1 The validity of equations ( 3 ), ( 4 ), and (5) can be demonstrated by 

referring to the limiting values of complex integrals taken along a path 
a short distance above the singular point on the real axis. In the case 
of equations ( 3 ) and ( 4 ) this process yields the Cauchy principal value. 
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In mathematical terms the problem is to minimize the double integral, 
equation (5), while holding fixed values of 


p+S 

L = pV T dy 

s 


and 


B = pV 



(y-s 0 ) r dy 


(O 


(7) 


where B is the bending moment specified about the point s D . For the 
time being s 0 will be taken as the origin, or wing root (s 0 = 0), 
although later another example will appear. 

Although the derivation of the formulas for induced drag makes use of 
the concept of the lifting line, it is important to note that the results 
are not actually restricted to this approximation. According to Munk*s 
well-known stagger theorem the induced drag of a lifting surface will be 
equal to that of a lifting line if the spanwise load distributions are 
the same. 

It should be noted further that the induced drag of a wing having a 
given lift and a given spanwise load distribution is not affected by the 
compressibility of the air at subsonic speeds. At supersonic speeds an 
additional drag associated with the formation of waves arises and the 
induced drag, which is associated with the vortex wake, becomes only a 
part, of the total pressure drag. 


THE DISTRIBUTION OF DOWNWASH FOR MINIMUM DRAG 


In general, if the drag is to be a minimum, a small variation in the 
shape of the curve of spanwise loading will produce no first-order change 
in the drag. The variation in shape may take the form of a small addition 
to the original loading; it is then necessary to find conditions under 
which the drag added by a small additional loading is zero. 

The solution of this latter problem is rendered especially simple 
by the mutual drag theorem (reference l), which arises from the evident 
symmetry of the integral to be minimized (equation (5))* The theorem 
states that if the lift distribution (represented by T (y)) is the sum 
of two distributions Ti, and r 2 , the drag of Tx arising from the 
downwash field of r 2 is exactly equal to the drag of arising from 

the downwash of r 1# 
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Consider now an initial distribution designed to achieve minimum 
drag. (See fig. 1.) The drag added by a small additional loading will 
be composed of three parts, namely: 

1. The drag of the additional lift acting alone 

2. The drag of the original loading arising from the dcwnwash 

field of the additional loading 

3. The drag of the additional loading induced by the downwash field 

of the original loading 


Item 1 is of second order in terms of the magnitude of the added lift for 
smooth distributions, that is, so-called "weak variations." (The fact 
that this second— order term is invariably positive insures that the drag 
will be a minimum and not a maximum. ) Items 2 and 3 are equal by the 
mutual drag theorem. The first-order variation in drag can then be com- 
puted by considering only the drag of the small additional lift acting in 
the induced downwash field wj_ ( y ) of the original lift. 

The conditions of fixed bending moment and fixed total lift are met 
by allowing only those curves of lift variation that produce no change in 
these quantities, that is, curves having zero area and zero moment. It 
can be seen that such curves of variation must have at least three ele- 
ments to meet tl^e conditions of zero area and zero moment. Furthermore, 
any curve meeting these conditions can be subdivided into groups of three 
elements so that the individual groups also satisfy the conditions. 

Hence, as the representative of such restricted curves of variation we 
may adopt three small elements having areas Z 2 , and l 3 (fig. 1). 

These elements, together with their positions y 2 , and y 3 and the 

local values of the downwash Wi x , etc., due to the original loading 
must satisfy the following three equations: 


for 

5L = 0, 

l l 4* ^2 ^ 3 = 0 

for 

5B = 0, 

l iJi + T-Z72 + l 3 y 3 = o 

for 

5Di = 0, 

1 l*ii + Z 2 W i 2 + 


It can be seen that these equations will be consistent if wj_ x ~ a+byi, 
wi 2 ~ a+by 2 and wj. 3 ~ a+by 3 , where a and b are constants to be 
determined from the given conditions. Since such equations must be. satis- 
fied for all positions yi, y 2 , etc., it is concluded that, in general. 


wi ~ a + by 


(9) 



5 


NACA TN 224-9 


Hence, for a minimum induced drag with a given total lift and a given 
bending moment the downvash must show a linear variation along the span. 2 
(See fig. 2. ) 

The foregoing method may be readily extended to a more general class 
of problems involving bending moments or rolling moments. Suppose, for 
example, a braced wing is considered, as in the dotted outline of fig- 
ure 3* In this case the bending moment developed by that portion of the 
lift acting inboard of the point of bracing attachment may be of no con- 
cern, but it may be desired to limit the bending moment developed by that 
portion of the spanwise load curve extending between this point and the 
tip. In this case s Q will not be zero. At least three elements are 
required to preserve stationary values of the lift and bending moment, 
and it is evident that at least two of the elements must lie to the right 
of the point s Q . The three simultaneous equations are (see fig. 2): 


ii + Z 2 + I 3 = 0 

Z2(y2-s 0 ) + ^3(^3“- s o) “ 0 

l lVp + Z 2 w-j + Z3W-1 = 0 

1 s X 3 


^ ( 10 ) 
4 


Here y 2 and y 3 are to the right of the point s D and yi lies to the 
left of this point. For these equations 3 to be consistent wi must have 
the form 


v il ^ ~ a+b(y 2 -s 0 ) ; vi 3 ^ a+b(y 3 -s 0 ) 


Hence , in general , the downvash will be a constant over the portion of the 
span for which the moment is not specified, as illustrated in figure 3- 
If no restriction whatever is placed on the moment there is obtained the 
solution of Munk*s original problem, namely, that the downvash should be 
constant over the entire span. 


It may be noticed at this point that, whereas the discussion has empha- 
sized the idea of minimizing the drag, the analysis actually makes no 
distinction between the lift, bend i ng moment, or drag, in that station- 
ary values of all three are demanded. Thus equation (9) may be consid- 
ered a necessary condition for the solution of the following problems: 
(l) given the total lift and the induced drag to find the distribution 
of lift over the span that will result in a minimum bending moment, 

(2) given the bending moment and the induced drag to find the dis- 
tribution resulting in the maximum total lift. 

3 See reference 3 for a discussion of solutions of such equations. 
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Determination of Span Loading and Induced 
Drag From the Down-wash. Distribution 


The case of bilateral symmetry with moment specified about the 
root section will serve as an example of the calculation of the actual 
span loading and induced drag. It will be evident from the foregoing 
that the downwash distribution will consist of two straight-line 
segments with a reversal of slope at the plane of symmetry. It is then 
necessary to compute the spanwise variation of f corresponding to 
such a curve of downwash. 


To perform this calculation by standard methods of airfoil theory, 
use is made of the idea that at a great distance behind the wing the 
vortex sheet forms a two-dimensional field of motion, with the discon- 
tinuity in the lateral velocity across the sheet given by dr/dy, and 
the downwash w given by twice the value of the induced downwash wp 
at the wing. Hence, the quantity l/2 (dr/dy) — 2iwp can be evaluated 
by means of the familiar complex velocity function v — iw of the two- 
dimensional potential theory using for v its value just above the 
vortex sheet. In this theory if the vertical component of velocity w 
is given along the line representing the trace of the span, then the 
velocity vector at any other point in the field § = y + iz may be 
obtained from the relation (reference b) 


= I 1 p+s ^(l) a/~s 2 — TjS 

’ V s 2 - ? 2 j_g i-e 


dq 


As noted above. 


dr 


^ = v(y + oi) — v(y — oi) = 2v(y + oi) 


( 11 ) 


( 12 ) 


Introducing w = a + by for y > 0 and w = a - by for y < 0 into 
equation (ll) yields, after integration. 


— = -2a 

ay 


/ 5 = 


iiq 

— b 
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and hence 

r = 2 (a + y S £ -y 2 + y 2 cosh -1 ~ 

The spanwise loading thus contains the elliptical distribution as one 
component. . - 
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Equation (13) for the spanwise distribution of circulation enables 
the determination of the over— all lift, bending moment, and drag in terms 
of the unassigned constants a and b. The use of equations (2), (6), 
and (7), together with the wing semispan s, yields the following values: 



► 




(HO 


It is convenient to specify the bending moment of the lift in terms of 
the lateral position of the centroid, or center of pressure, of the load 
curve. The lateral centroid as a fraction of the semispan s may be 
denoted by y 1 (i.e., y 1 = 2B/Ls). Then, solving for a and b. 



The expression for induced drag In terms of the lift and the lateral 
center of pressure becomes 


D i 


n | V 2 (2s) 2 



rt 2 y* 2 — 12rty' 



( 16 ) 


This equation yields the minimum drag for the given position of y' . If 
the lateral center of pressure is specified so as to coincide with that 
for an elliptical loading (i.e., b =0; y’ = ^/3 rt )> then the above 
formula reduces to 



L 2 

jr | V 2 (2s) 2 


(17) 


The optimum distribution of loading for a given position of the centroid 
y T may be obtained from equation (13) with the aid of equations (15). 
The result is 
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Drag For a Given Bending Moment With Unrestricted Span 


The foregoing calculations show, as was to he expected, that the 
elliptic loading yields a smaller drag than any of the others within a 
restricted span. However, if the restriction on the span is removed, 
still lower values of the induced drag can he obtained without any 
increase in the bending moment at the wing root. The lower values are 
obtained by permitting the span to increase and at the same time adopt- 
ing a more tapered form of the loading curve. 


Equation (1 6) which contains the three variables lift, span, and 
center of pressure can be easily rearranged to show the variation of 
drag with span when the bending moment and the lift are held at fixed 
values. In this case, the lateral position of the center of pressure 
y T s will be fixed, while the form and extent s of the load curve will 
vary. In order to provide a convenient basis for comparison the span 
and shape of the load curves will be related to the elliptic loading. 

If s/s e denotes the ratio of the semispan of the wing to that of an 
elliptically loaded wing having the same total lift and bending moment, 
then equation (1 6) can be rewritten: 


D, = 


* | V 2 (2s e ) 2 


*©‘ -“(»'*<»■] 


(19) 


The quantity in the bracket is the ratio of the induced drag to that of 
the corresponding elliptically loaded wing. This ratio is plotted in 
figure k to show the decrease of drag possible by increase of the span. 
The forms of load curve required for the minimum drag at various values 
of s/s e are shown In figure 5* 


It will be noted that a 15—percent reduction of the induced drag 
below that for elliptic loading can be achieved with a 15— percent 
increase In span. Further increases of span between 15 percent and 
50 percent (s =1.15 to 1.50) yield no significant reductions, however. 

At still larger values of s the drag becomes lower, and approaches zero 
at an infinite value of s. For extreme values of s/s e the curves begin 
to show negative loadings at the tips and eventually the bending moment 
at certain points along the span will exceed that at the wing root. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Sept. 25, 1950. 
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APPENDIX 

DEFINITIONS OF SYMBOIS 

L total lift 

l element of lift 

induced drag 
B bending moment 

p air density 

r circulation 

w^ induced downwash velocity at wing 

w downwash velocity, at infinity ( w = 2w^) 

v lateral velocity 

V velocity of flight 

y,Tj distances along wing semispan 

s c point of origin for bending moment 

s length of wing semi span 

y* lateral position of load centroid as a fraction of s 
a,b constants 
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FIGURE I- SPANWISE LOAD CURVE WITH THREE 
ELEMENTS OF VARIATION 



FIGURE 2- DOWNWASH DISTRIBUTION FOR 
MINIMUM DRAG WITH RESTRICTED BENDING 
MOMENT ABOUT CENTER SECTION 
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FIGURE 4.- VARIATION OF INDUCED DRAG WITH SEMISPAN RATIO 
S/S e , FOR WINGS HAVING A FIXED TOTAL LIFT AND A 
FIXED BENDING MOMENT. 
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FIGURE 5.- VARIATION OF SHAPE OF THE SPANWISE LOADING 
CURVE WITH SEMISPAN RATIO S/S e , FOR WINGS HAVING 
A FIXED TOTAL LIFT AND A FIXED BENDING MOMENT 




